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Abstract 

We exploit a recent lattice investigation (UKQCD) on the topological structure of the 
(quenched) QCD vacuum, in order to gain information on crucial building blocks of instanton 
perturbation theory. A central motivation is to further constrain our previous predictions 
of instanton-induced hard scattering processes. First, we address the generic problem of 
extracting quantitative information from cooled lattice data. We find a new scaling variable, 
interpreted as a "cooling radius", which allows to combine lattice data for a whole range 
of lattice spacings and cooling sweeps. This variable strongly helps to extract information 
on the uncooled distributions of interest. After performing the continuum extrapolation 
of the instanton size and instanton-anti-instanton distance distributions, we find striking 
agreement with the theoretical predictions from instanton-perturbation theory, for instanton 
sizes <0.5 fm and distances >0.5 fm, respectively. These results imply first direct support 
for the validity of the known valley interaction between instantons and anti-instantons. 



1. Non-abelian gauge theories like QCD are known to exhibit a rich vacuum structure. The latter 
includes topologically non-trivial fluctuations of the gauge fields, carrying an integer topological 
charge Q. The simplest building blocks of topological structure are instantons (Q = +1) and 
anti-instantons (Q = —1) which are well-known explicit solutions of the euclidean field equations 
in four dimensions 0. 

Instantons are widely believed to play an important role in various long-distance aspects |2| of 
QCD: 

First of all, they may provide a solution of the famous Ua(1) problem || (m v > ^> m v ), with the 
corresponding pseudoscalar mass splitting related to the topological susceptibility in the pure 
gauge theory by the well-known Witten-Veneziano formula ||J]. Moreover, a number of authors 
have attributed a strong connection of instantons with chiral symmetry breaking as well as 
the hadron and glueball spectrum. 

However, there are also very important short- distance implications [§]HTIJ of QCD instantons: 



They are known to induce certain processes which violate chirality in accord with the general axial- 
anomaly relation || and which are forbidden in conventional perturbation theory. Of particular 
interest in this context is the deep-inelastic scattering (DIS) regime. Here, hard instanton-induced 
processes may both be calculated []8| [Toll within instanton-perturbation theory and possibly be 



detected experimentally [|7yn|,l2"|,r3|. Since an experimental discovery of instanton-induced events 
would clearly be of basic significance, our ongoing theoretical and phenomenological study of the 
discovery potential of instanton-induced DIS events is well motivated. 

On the other hand, there has been much recent activity in the lattice community to "measure" 



topological fluctuations in lattice simulations |TJ| of QCD. Being independent of perturbation 



theory, such simulations provide "snapshots" of the QCD vacuum including all possible non- 
perturbative features like instantons. They may be exploited to provide crucial support for 
important prerequisites of calculations in DIS, like the validity of instanton-perturbation theory 
and the dilute instanton-gas approximation for small instantons of size p < 0(0.3) fm. Along 
these lines, we were able |9|] to translate the lattice constraints into a "fiducial" kinematical region 
for our predictions of the instanton-induced DIS cross-section based on instanton-perturbation 
theory. 

The purpose of the present letter is to present a more thorough confrontation of predictions from 
instanton-perturbation theory with recent lattice simulations. Specifically, we concentrate on 
recent high-statistics results from the UKQCD collaboration fl5|]16l , try to carefully discuss the 
effects of cooling and perform an extrapolation to the continuum before drawing our conclusions. 

2. Let us start by considering various important quantities relevant for instanton (J)-induced 
scattering processes, that may both be calculated in /-perturbation theory and measured in 
lattice simulations. In /-perturbation theory one expands the path integral for the generating 
functional of the Green's functions about the known, classical instanton solution = A^p + . . . 
instead of the trivial vacuum A£> = (pQCD) and obtains a corresponding set of modified 
Feynman rules. Like in conventional pQCD, the strong gauge coupling a has to be small. 

The /-induced contribution to hard scattering processes is typically encoded in the total cross- 
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section fl,|lO|l associated with an /-induced partonic subprocess q' p =^> X, 



OO OO / 2 — \ 

a {I \q'p^X) = Jd 4 R JdpJdpD(p)D(p)JdUe~^ n V^ ,Ep ){...} 



o o 



a'- 7 '' involves integrations over all / [7]- "collective coordinates", including the / [J]-sizes p [p], the 
JJ-distance[| 4- vector and the II relative colour orientation U . With each initial-state parton, 
there is an associated "form factor" MM . 



^7r/ (2 x) exp(— x) for x — > +oo; 
/'(.'•) = .rh\(.r) { 1 for x = 0; (2) 

oscillating for x imaginary. 

An important quantity entering Eq. ([!]) as generic weight is the (reduced) /-size distribution in 
the dilute gas approximation, D(p). It is known in the framework of /-perturbation theory for 
small a(p r ) ln(p/z r ), where p r is the renormalization scale. After its pioneering evaluation at 
1-loop |J for N c = 2 and its generalization ]T7| to arbitrary N c , it is meanwhile available |TB| in 
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2-loop renormalization-group (RG) invariant form, i.e. D 1 dD/dln(p r ) = 0(a 

D( P ) = %(^) 2 ^ xp (-^) (p , r) «. + <*-".*>^, (3 ) 

^ - ^Sv^ e "' 5n371 " +a29 ' 746 "' (Ref ' |19|)i (4) 

For the case n./ = 0, relevant for our comparison with lattice results below, the reduced size 
distribution D(p), Eq. @, equals the proper size distribution, i.e. the number of instantons per 
unit volume per unit size, dni/d 4 x dp, (in the dilute gas approximation). On the other hand, for 
rif > 0, the proper size distribution reads 

~ D m {p) = D(p) lliprrnM) (p^' 70 ^^, (6) 

for small prrii(p r ), where rrii(p r ) are the running quark masses. However, even for n/ > 0, only 
the reduced size distribution D(p) appears in cross-sections like Eq. ([]]), since for small rrii, the 
mass-dependent factor in Eq. (|6]) is cancelled by corresponding terms from the external quarks. 

^oth an instanton and an anti-instanton enter here, since cross-sections result from taking the modulus squared 
of an amplitude in the single /-background. Alternatively, one may view the cross-section (Q) as a discontinuity 
of the q'p forward elastic scattering amplitude in the //-background || . 
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The powerlaw behaviour of the (reduced) J-size distribution, 

D(p)~/°- 5 +°W, (7) 

generically causes the dominant contributions to the J-size integrals to originate from the in- 
frared (IR) regime (large p) and thus often spoils the applicability of /-perturbation theory. In 
deep-inelastic scattering, however, one parton q' say, carries a spacelike virtuality, such that the 
contribution of large instantons to the integrals ([!]) is exponentially suppressed (c. f. Eq. (0)), 



F(^-q' 2 p)F^-q' 2 p) (xe- Q ' {p+p) - for Q' = ^ -q' 2 > 0, (8) 

and the predictivity of /-perturbation theory is retained for sufficiently large Q'. 

A second important ingredient into Eq. (TI]) is the function fl(U,R 2 /(pp),p/p), appearing in the 
exponent with a large numerical coefficient 4ir/a. It incorporates the effects of final-state gluons. 
Within strict /-perturbation theory, it is given in form of a perturbative expansion P,p0f for 
large R 2 , while in the so-called //-valley approximation^ |]2T| is associated with an analytically 
known closed expression^] |^3| for the interaction between / and I, Q— a/(4:7r)S[A 11 ] — 1, 

. .), f , A = (± - 2 ,) (1 + , 2 ) E X + (Q - 2 Z + 4 *>) (1 + 2 , 2 ) S 2 (9) 




Si 



l-A 4 +4A 2 ln(A) . y l-A 2 + (l+A 2 )ln(A) ^ 

(A 2 -l) 3 ' 2 (A 2 -l) a ' 



(10) 



pp p p 



The //-valley interaction (|9|) only depends on two of the four angle variables necessary in general 
to parametrize the relative color orientation U (see e.g. Ref. f24j). In terms of these variabels 
(v = cos a, z = cos(5) 2 ) and an additional azimuthal angle variable 0, the integration over the 
SU (3) group measure then takes the form |TD 



du=- [ \ V \ (i-v 4 )dv [ d(f>r 2 ^ — %- — , = i. (ii) 

7T 7-1 ' ' V J J Y J 1 _ 2^COS(0) + V 2 V ' 

With these prerequisites in mind, let us next consider the group-averaged distribution of // 
pairs, drij-j/ d i x d A Rdp d~p. This quantity is of basic interest for /-induced scattering processes 
([]]), crucially depends on the //-interaction Q and may be measured on the lattice. For the pure 
gauge theory^] (n/ = 0) and in the dilute-gas approximation, it takes the form |£,10] 



drijj 



d A xd A Rdpdp 



D l7 (R, p, p) = D{p) D{p) J dU exp 



a- 



4 * -afuXX 



ms( s /j/ Vf>P) \ ' PP' P, 



(12) 



2 For the most attractive relative //-orientation the form of the //-valley action was first given in Ref. |2^] . 
3 For rif > 0, there is an additional fermionic interaction factor | J d A x k oi (x) [i fi^ 11 ^] k j(x — R) \ 2nf under 



the group integral in Eq. (112), with kq being the fermionic zero modes. 
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with the scale factor Sjj = 0(1) parametrizing the residual scheme dependence. 

In the deep-inelastic regime, relevant for instanton searches at HERA, the collective coordinate 
integrals involved in the hard cross-section ([I]) are known to be dominated by a unique saddle 
point [^]. Therefore, the Bjorken variables (Q',x') of the /-induced subprocess in e ± P scattering 
are effectively related in a one-to-one correspondence to the collective coordinates as (c. f. also 
Eqs. ©, (§)), 

p ~ p ~ l/Q' and R 2 ~ l/(p + q') 2 pp/R 2 ~ (p + </) 2 /Q' 2 = - 1. (13) 

Restrictions on the collective coordinates obtained from confronting predictions of /-perturbation 
theory with lattice results may thus be directly translated into a fiducial kinematical region for the 
predicted rate of /-induced processes at HERA 0. This fact illustrates the great importance of 
lattice simulations for /-induced scattering processes and largely motivates the present analysis. 

3. Before turning to our concrete analysis, we first have to address various topical problems 
associated with lattice data on instantons. The recent quenched QCD results (n/ = 0) from the 
UKQCD collaboration [p^Jl^l , which we shall use here, involve a lattice spacing a = 0.05 — 0.1 
fm and (roughly constant) volumes V = Z 3 • Z t ime of 32 3 64 a 4 at /3(a) = Q/ gl(a) = 6.4, 24 3 48a 4 
at $ = 6.2 and 16 3 48 a 4 at (5 = 6.0. 

In principle, such a lattice allows to study the properties of an ensemble of topological gauge 
field fluctuations (/'s and Z's) with sizes a < p < V 1 ^. However, in the raw data, ultraviolet 
fluctuations of wavelength 0(a) dominate and completely mask the topological effects which 
are expected to be of much larger size. In order to make the latter visible, a certain "cooling" 
procedure has to be applied first. Cooling is a technique for removing the high-frequency non- 
topological excitations of the gauge-field, while affecting the topological fluctuations of much 
longer wavelength /)>a comparatively little. After cooling, Z's and Fs can clearly be identified 
as bumps in the topological charge density and the Lagrange density. By applying sophisticated 
pattern-recognition algorithms, various distributions characterizing the ensemble of Fs and Fs 
may be extracted. 

However, during cooling instantons are also removed; either if they are very small (lattice artifacts) 



or through // annihilation. In Ref. |i~5 |, so called under- relaxed (or slow) cooling is used to reduce 
the latter problem. Generically, the cooling procedure introduces substantial uncertainties in the 
extraction of physical information on the topological structure of the QCD vacuum. No objective 
criterium is yet available on the quantitative amount of "vacuum damage" introduced by it (see 
Refs. 00). 



Since one is interested in the physics of the uncooled vacuum, it appears suggestive (at first) to 
try and extrapolate [T3J] the measured distributions to the limit of a vanishing number of cooling 



sweeps, n coo i — > 0. However, to do so would be to ignore the fact that the pattern-recognition 
procedures become increasingly unreliable near that limit, since the extracted instanton density 
appears to strongly increase for decreasing ra coo i. Below a certain minimal number of cools, 
(^cooi)min , even the total topological charge Q = nj — rij starts to deviate from the plateau value 
it has for a whole range of larger n coo i- 
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In view of these objections, we follow the alternative strategy of Ref. [15] to extract information 
about the physical vacuum. The method incorporates certain scaling studies that are important 
for the continuum limit a — > and consists in the following steps: First, at /3 — 6.0, say, one picks 
the smallest number of cooling sweeps for which the pattern-recognition starts to become reliable 
and the cooling history of the total topological charge enters a plateau. For straight cooling, one 
finds (n coo i) min ~ 10 corresponding to ~ 23 under-relaxed cools fllBl . The next step is to search 
for further, so-called equivalent pairs (f3, n coo {), for which shape and normalization of specific 



distributions (expressed in physical unitsQ) essentially remain the same (see also Ref. f25|). The 
result of Ref. [|15|] were the pairs (/3,n coo {) = (6.0,23), (6.2,46), (6.4,80). With these equivalent, 
scaling data sets, the continuum limit is then performed with hopefully little "vacuum damage" 
due to cooling. 

Yet, the crucial question remains at this point, whether the choice (/3,n coo {) = (6.0,23) as a 
"starting pair" for the equivalent data is somehow distinguished, i. e. physically meaningful. In 
this context, let us present an intriguing line of arguments: 

"Equivalent" data should first of all be characterized by a similar amount of cooling. For n coo i 
cooling sweeps, an effective "cooling radius" r coo i in physical unitsB may be defined, via a (quali- 
tative) random walk argument, as 



'"cool 

2^ 



const {^/n 



cool J 



a+5 



2r 



5 small. 



(14) 



The three equivalent pairs (/?, n coo \) above, indeed, turn out to have a very similar value of r coo i 
(see Table [I]). 





(6.0,23) 


(6.2,46) 


(6.4,80) 


y/n coo \ ^ 


0.447(2) 


0.460(2) 


0.459(2) 



Table 1: The "cooling radius" for the equivalent data sets is similar. 



Amazingly, the usefulness of the variable ( |I~4"D seems to be much greater. In Fig. [I], we illustrate 
that the average instanton size (p) in physical units from Ref. [fTq] , for all available values of 
(/3, n coo i), fits onto a single smooth curve as function of r coo i in Eq. ( |I~4"D with 6 small. In contrast, 
one obtains three widely spaced curves, if (p)/2r is displayed as a function of n coo i alone for fixed 
(3 = 6.0, 6.2 and 6.4, respectively. The same scaling with r coo i holds practically for all quantities 



extracted in Ref. ||15|| , characterizing the instanton ensemble on the lattice. 

Perhaps the most striking feature of Fig. [I] is the clear indication of a plateau corresponding to 
(p) ~ 0.5 fm. Let us try to provide some interpretation for it. Suppose, the localized fluctuations 
on the lattice are a mixture of two distinct ensembles, separated by a certain gap in their respective 
average sizes. For example, one expects fluctuations that are lattice artifacts typically located 



4 For consistency and minimization of uncertainties, one should use only a single dimensionful quantity to relate 
lattice units and physical units. Throughout our analysis, all dimensions are therefore expressed by the so-called 
Sommer scale j2q,p7| tq, with 2ro ~ 1 fm, which we prefer over the string tension fll5f. 
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Figure 1: The average instanton size from Ref. [15] in physical units for (3 = 6.0 [□], 6.2 [o], 6.4 [A] 
and for all available values of n coo i as a function of the effective cooling radius (0). There is a 
clear indication of a plateau corresponding to (p) ~ 0.5 fm, which is interpreted in the text. The 
quantity r denotes the Sommer scalei with 2 r — 1 fm. 



around smaller sizes (9 (few -a), while the physical instanton sizes are peaking around p ~ ~ 0.5 
fm. By gradually increasing the cooling radius r coo i, we first filter out the smaller artifacts, leading 
clearly to a corresponding increase in the effective (p). When the radius of our filter further 
increases, one then expects a plateau in (p) to indicate that the small artifacts have been mostly 
wiped out, while the much larger physical instantons are not yet strongly affected. Eventually, 
the instantons will also start to be erased, leading again to an increase in (p). 

In conclusion, we think that the plateau region is the correct regime for extracting physical 
information about instantons. The equivalent data sets of Ref. |15| are precisely located around 
the plateau. 

Our scaling variable r coo i allows to combine lattice data for a whole range of lattice spacings a 
and cools n coo i, which strongly helps to increase the separation power between the a-dependent 
artifacts and the instantons of a certain physical size. In future analyses, it may be worthwhile to 
determine the unknown constant in the definition ( |l~4T) of r coo i, e.g. by means of calibration via 
cooling an isolated instanton until it disappears. 

4. We are now ready to discuss the continuum limit of the UKQCD data for the size distribution 
of I's and J's, dn I+ j/d A x dp, and to compare it with the dilute-gas prediction D I+ j(p) = 2D(p) 
of /-perturbation theory, Eq. ([|), for rif — 0. 

We note, first of all, that the (small) residual a dependence of the equivalent data for the size 
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distribution in physical units may be parametrized as 



^M =baaSm {jh)- (15) 

The continuum limit of Eq. (|T5f) can thus be performed quite reliably, by simply rescaling the 
arguments p — > (p(0))/(p(et)) • p. Here, (p(0))/2r = 0.518(5) denotes the continuum limit of the 
weakly varying average p values, (p(a)). As in Ref. [115] , (p(a))/2r was extrapolated linearly in 
(a/ro) 2 . The resulting continuum size distribution obtained by the above rescaling from the three 
equivalent data sets is displayed in Fig. ||. It scales nicely. 
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Figure 2: Continuum limit of "equivalent" UKQCD data [119 for the (J + J)-size distribution 
at (/J.nooob) = (6.0,23) [□], (6.2,46) [o], (6.4,80) [A]. The striking agreement with 2D(p) of 7- 
perturbation theory from @ is apparent. The 3-loop form of aj^ with Aj^g n =0 from ALPHA f28 | 
was used, (a) Log-log plot to exhibit the expected power law ~ p 6 and the agreement in magnitude 
for small p over a wide range of p r . The dashed error band results from varying A^g n ^ =0 and p r 



within its error and given range, respectively; (b) For p r 
the peak. 



1.2/p, the agreement extends up to 



We are now ready to perform a quantitative comparison with /-perturbation theory. The corre- 
sponding prediction @ is a power law for small p, i.e. approximately D ~ p 6 for rif = 0. Due 
to its 2- loop RG-invariance the normalization of D I+ j(p) is practically independent of the renor- 
malization scale p r over a wide range. It is strongly dependent on tq Aj^g, for which we take the 
most recent, accurate lattice result by the ALPHA collaboration |28[] , 2r A^g n ^ =0 = (238 ± 19) 
MeV fm. In Fig. Q (a) we display both this parameter-free prediction (^) of /-perturbation theory 
and the continuum limit of the UKQCD data in a log- log plot, to clearly exhibit the expected 
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power law in p. The agreement in shapef] and normalization for p/2r <0.3 is striking, indeed, 
notably in view of the often criticized cooling procedure and the strong sensitivity to Aj^g. One 
may argue that a choice pp r = 0(1) is theoretically favoured, since it tends to suppress ln(pp r ) 
corrections. Indeed, as apparent from Fig. [2] (b), the agreement extends in this case up to the 
peak around p/2r ~ 0.5. 

At this point, an interesting way to proceed is to define a non-perturbative "instanton-scheme" 
for a by the requirement that the theoretical expression (|3|) for the size distribution D I+ j(p) 
precisely reproduces the corresponding lattice data, 

By comparing in the perturbative regime with the MS form (H) for p r = s/p, we arrive at a 
standard scheme conversion formula, 



with 



«ms(*/p) = «i(Vp) - sz Ms/ si) «i(l/p) 2 + 0(«i(l/p) 3 ), (17) 



Sl = ^MS and a m (s 1 /p)~a 1 (l/p). (18) 
Al 

Next, we note that as a function of aj, D I+ j(p) has a peak structure like the lattice data, with 
maximum at 

7T 2 

iVc l+v /l + i ln ^i)(2/?o-^) 

The scale factor si of the /-scheme is non-perturbatively and uniquely determined by the peak 
value of the lattice data for the size distribution. Using the very good Gaussian fit to all data in 
Fig- 1(b), 



5 dn I+1 (qp) 6 exp(-(cp) 2 ) . b = 1.683(15)/2r 

1 0) d*xdp l&t (l + (bp)P) ' c = 2.525(52)/2r 



4.908(40)/2r 
1.683(15)/ 
2.525(52) 7 
12.43(25) 



we find 



si = 1.18 and a I]max = 0.98. (21) 

In the perturbative regime, the /-scheme is thus quite close to the MS-scheme (c. f. Eq. (|T8"D). 

We are now able to extract for each data point of the / + / size distribution a corresponding value 
of ai, as illustrated in Fig. || Actually, besides the asymptotically free solution [ai] to Eq. (PH), 
there is a second one [a[ 2 ' ) ] (dotted line in Fig. §), with a monopole-like relation to oti, 

"i (2) ■ ai = 0(1), (22) 

5 The agreement in shape up to p ~ 0.3 fm has been noted by us already previously j^]. Here, we have extended 
the analysis to the continuum limit and, in particular, verified the agreement in normalization. 
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Figure 3: (a) Running coupling ai in the "instanton" scheme, extracted via Eq. (^) from 
the I + I size distribution. The same UKQCD data as in Fig. |2] have been used. The dotted 
line denotes a possible second solution vanishing in the IR limit. The excellent agreement with 
the corresponding 3-loop perturbative a^(si/ p) ~ a\{l/ p) down to 1/p ~ 0.4 GeV is apparent 
(solid line), (b) At large distances, a\ from (a) appears to be quite close to a qq defined via the 
static qq potential, with (4/3r 2 )a, 
2r 0y /a = (472.4 ±4.3) MeVfm. 



qq 



1/r) approaching the string tension a for large r and p7 



vanishing in the IR-limit and increasing like af\p) oc (p/kY° /2Nc in the UV-limit. This obser- 
vation may be of interest in the context of Ref . |29[ . 

It is remarkable that in the non-perturbative long-distance regime, p/2r >0.5, a\{l/ p) from the 
I + I size distribution tends to approach a qq (l/p) defined via the static qq potential, Vq q , 

4 cx qq (l/r) _ dV m {r) it 1 

3 ar 12 r z 

with string tension a. This is illustrated in Fig. [I] (b). 

In summary, a\ as non-perturbatively defined through Eq. ([RJ), is close to the MS-scheme in the 
perturbative short-distance region and to the gg-scheme in the non-perturbative long-distance 
regime. Figs. [| (b) and |3] (b) illustrate that the transition between perturbative and confinement- 
sensitive physics seems to be rather sharp and located where the instanton density peaks, around 
p ~ 0.5 fm. 



The above results also suggest the following working hypothesis which might considerably 
enlarge the region of applicability of predictions like ([I]) and ( |T2| ) for observables associated with 
sufficiently dilute instantons. Whenever one encounters the /-size distribution D(p) within the 
perturbative calculus, one may try to replace the indefinitely growing perturbative expression 
-D(p) ~ pA>~ 5 for the /-size distribution by the functional form actually extracted from the lattice 
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(see e.g. Eq. (20)). Obviously, such an attempt eliminitates the problematic IR divergencies of 
/-size integrals. Due to the strong peaking of the measured size distribution (|20|), only J's of size 
p ~ 0.5 fm effectively contribute to the observables of interest^] and the dilute gas approximation 
may well continue to hold qualitatively up to the peak around 0.5 fm. 

5. As an immediate first application of this strategy, let us turn to the group-averaged distribu- 
tion of //-pairs, Djj(R, p,p), in the dilute gas approximation, Eq. ([121) . While the conservative 
prediction of instanton perturbation theory is restricted to p, p<0.3 — 0.5 fm, UKQCD lattice 
results jnjnjl are available up to now only for the //-distance distribution, integrated over all 
values of p,p, 

d 4 xd 4 R ~ J J P 9 d 4 xd*Rdpdp ' ( ' 

This is the point to exploit the above working hypothesis: We now replace in the theoretical 
expression for the //-distance distribution (|24f) along with Eq. (|12"D the perturbative form for 
D(p) by the Gaussian fit (pD|) to the UKQCD lattice data. We are left with only one free 
parameter, the scale factor Sjj = 0(1) in the running coupling a^(s T j / \fftp) . On account of the 
peaking of the size distributions in Eq. flT2]), the running coupling at a momentum scale around 
(9(0.3) GeV enters significantly. In this region, we make use of the non-perturbatively extracted 
coupling in the /-scheme (Fig. ^), with 



a m -£= )^aj M^-= ; sj = 1.18, (25) 




as a smooth continuation of the perturbative MS coupling (c. f. Eq. (0)). In the quantitative 
evaluation of Eqs. ([24]) and (0), we perform two integrals of the group integration and the two 
integrals over the /- and /-sizes numerically. 

Before being able to compare this theoretical prediction of the //-distance distribution with lattice 



results, we have to perform the continuum limit of the respective UKQCD data ||16|| . This we 
do in complete analogy to the method described in Sect. 4., Eq. (p~5|) , for the /-size distribution. 
Unfortunately, the equivalent data for j3 = 6.4 are missing in this case [ 16| . The equivalent data 



for the remaining two (3 values, ((3,n coo \) = (6.0,23) and (6.2,46), again are well compatible with 
a scaling law, like Eq. (|15|), 

" 1 n - function — - . (26) 



d 4 xd 4 R I d 4 xd 4 R\ R large V r ( a ), 

By extrapolating the scaling factor A(a) = r(a)/r\p = Q 2 linearly in (a/ro) 2 , we find from requiring 
best overall scaling (p6|) of the two equivalent data sets, A 7 j(a) = 1.32 — 69.75(a/2r ) 2 . The result- 
ing "continuum data" obtained from the simple rescaling of the arguments, R r(0)/r(a) R = 
A(0)/A(a) R, are displayed in Fig. f| and exhibit good scaling. 

The solid line in Fig. [| denotes the theoretical prediction discussed above with the scale factor 
Sjj = 0.82 being 0(1) as expected. 

6 This philosophy is partly also inherent in the /-liquid model 0. 
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Figure 4: Continuum limit of "equivalent" UKQCD data for the //-distance distribution, 
drijj/cPx d 4 R, normalized to its value at R/2r = 1.5, at (/3, n coo i s ) = (6.0, 23) [□], (6.2, 46) [o]. As 
discussed in the text, the solid line denotes the theoretical prediction from instanton-perturbation 
theory based on the valley form of the //-interaction. 



We note a very good agreement with the lattice data down to //-distances R/(p) — 1, corre- 
sponding to £>3 (c. f. Eq. (|10D ). These results imply first direct support for the validity of the 
valley interaction (|9|) between instantons and anti-instantons. 
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